Entanglement and bifurcations in Jahn- Teller models 



^1" 
O 

o 

(N 

X> 

<D 
IX, 

(N 
> 

o 

<N 
O 

o 

9 L,: 
^— > ■ 

c ■ 

=3 : 
cr 



X 



Andrew P. Hines0 Christopher M. Dawson, Ross H. McKenzie, and G.J. Milburn 

Centre for Quantum Computer Technology, School of Physical Sciences, 
The University of Queensland, St Lucia, QLD 4072, Australia 
(Dated: February 9, 2008) 

We compare and contrast the entanglement in the ground state of two Jahn- Teller models. The 
E ® f3 system models the coupling of a two- level electronic system, or qubit, to a single oscillator 
mode, while the E (g> g models the qubit coupled to two independent, degenerate oscillator modes. 
In the absence of a transverse magnetic field applied to the qubit, both systems exhibit a degenerate 
ground state. Whereas there always exists a completely separable ground state in the E®f3 system, 
the ground states of the E Cg> e model always exhibit entanglement. For the E ® f3 case we aim 
to clarify results from previous work, alluding to a link between the ground state entanglement 
characteristics and a bifurcation of a fixed point in the classical analogue. In the E®e case we make 
use of an ansatz for the ground state. We compare this ansatz to exact numerical calculations and 
use it to investigate how the entanglement is shared between the three system degrees of freedom. 



I. INTRODUCTION 

The burgeoning held of quantum information science 
has provided new tools with which to probe the charac- 
teristics of complex quantum many-body systems. More 
specifically, the study of the entanglement properties of 
systems is an active area of research, aimed at shedding 
new light on previously studied fundamental systems. 

As such there has been numerous studies of the en- 
tanglement in the ground states of various systems (see 
H H II II H H toll IS 113 and reference therein). Of 
particular interest has been those systems which exhibit 
a quantum phase transition, where it has been demon- 
strated that the entanglement properties are connected 
with this critical phenomena 0, 0, El El • 

Another problem where an understanding of the entan- 
glement properties offers a new perspective is in the study 
of decoherence. Any real-life quantum system interacts 
and becomes entangled with its environment, causing 
quantum superposition states to decohere into classical 
statistical mixtures. One way of studying the process of 
decoherence in open quantum systems is by the quan- 
tum environment and studying the now closed system- 
environment setup. 

Probably the most well-known system-environment 
model is the spin-boson model [T^. Il3|. which describes 
the interaction between a qubit (any two-level system) 
and an infinite collection of harmonic oscillators, mod- 
eling the environment. The entanglement between the 
qubit and its 'environment' (the oscillators) in the ground 
state of this model was recently studied by Costi and 
McKenzie Q where a further link between entanglement 
and QPT's was established. 

As a way of investigating the decoherence induced by 
certain measurements, Levine and Muthukumar |14| con- 
sider a model describing a qubit coupled now to a single 
environmental mode. This system is also known as the 
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E ® (3 Jahn- Teller model [l5| . Levine and Muthukumar 
|l4j study the variation in the ground state entanglement 
with respect to the strength of the coupling between the 
qubit and the oscillator. In the massive limit (m — > oo) of 
the oscillator, two parameter regions are identified, where 
the ground state is completely separable and where the 
qubit and oscillator are entangled. In this article we aim 
to clarify this result, in the light of previous results from 
the authors and Lambert, Emary, and Brandes |||, 
regarding ground state entanglement and corresponding 
fixed point bifurcations in the classical analogue. 

Following the natural progression from the single os- 
cillator case, we consider the E (g) e Jahn- Teller system, 
which describes the coupling of a qubit to two identi- 
cal (uncoupled) oscillators. Jahn- Teller models are of 
great importance in the study of the geometry of molec- 
ular structure, in cases where the coupling between elec- 
tronic and nuclear states cannot be ignored |30|. The 
E (g) e Jahn- Teller system describes the coupling between 
a doubly degenerate electronic state (E) and a doubly de- 
generate normal mode (e). Such a model has been used 
to study the degree of electron-nuclear entanglement in 
molecular states [r7^. 

In the case of both Jahn- Teller models considered here, 
when there is no transverse magnetics field applied to the 
qubit, the ground state has a two- fold degeneracy. This 
means that there are an infinite number of ground states, 
consisting of all possible superpositions of any two or- 
thogonal ground states. Not all such ground states will 
necessarily contain the same amount of entanglement. To 
obtain a complete picture of the ground state entangle- 
ment, one has to consider the entanglement in all possible 
ground states. 

Often, it is the case that there simultaneously exists 
ground states with maximal entanglement, and com- 
pletely separable ground states. Certainly, it can be 
shown that for a system of two qubits, if there are two 
orthogonal, maximally entangled ground states, then an 
equal superposition of the two is completely separable. 
This is the case for the E ® (3 model, where, irrespec- 
tive of the strength of the coupling, there is always a 
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ground state which contains no entanglement. However, 
the E®e model exhibits the intriguing property that for 
all ground states when the coupling is greater than zero, 
the qubit is entangled with the oscillators. While in the 
limit of large coupling there are ground states with maxi- 
mal qubit-oscillators entanglement, we show that the en- 
tanglement in all ground states is always bounded below 
by some non-zero value. 

We begin with the E®(} model, by analyzing the corre- 
sponding classical model before considering the entangle- 
ment in the ground state. This is followed by the same 
analysis for the E ® e model and a comparison of the 
ground state entanglement characteristics of the two. 



II. E ® [3 : A QUBIT COUPLED WITH A 
SINGLE OSCILLATOR MODE 

The E ® (3 is the mathematically simplest Jahn- Teller 
effect, and occurs where a doubly degenerate state (the 
qubit) becomes coupled by a single boson mode (the 
oscillator). The entanglement characteristics of such a 
model system have been recently studied by Levine and 
Muthukumar Q , where they considered a qubit coupled 
to a single harmonic oscillator described by the Hamilto- 



H = Aa x + L 



[a + a> a. 



(1) 



where u> is the natural frequency of the oscillator, L is 
the coupling strength and A is the strength of the trans- 
verse magnetic field acting perpendicular to the coupling, 
all of which are in units such that h — 1 (for the rest of the 
paper we assume unit mass, m = 1). This Hamiltonian 
can also be written in terms of the position coordinate, 
q of the oscillator, as 



H = A<7 T 



Lqdr z 



if- 



2 -2 
Lu If 



(2) 



This system is a simpler version of that studied by 
Emary and Brandes 0], who considered a collection of 
N two-level atoms, modeled as a single collective spin, 
interacting with a single bosonic mode via a dipole inter- 
action - the so-called Dicke Hamiltonian. In their analysis 
based on functional integrals, Levine and Muthukumar 
|14| identified a critical parameter value corresponding to 
a qualitative change in the ground state of the system. In 
the next section we consider the analogue classical system 
and derive this critical parameter via a simple analysis 
of the dynamical fixed points. 



tion of the fixed points [l^ in the corresponding classical 
system. 

Letting q and p be the classical position and momen- 
tum coordinates of the oscillator, and L x ,L y and L z the 
spin coordinates of the spinning top (the classical ana- 
logue of the qubit), the equations of motion are found to 
be 

q = p, (3a) 
P = -LL Z - uj 2 q, (3b) 
L x = -LqL yi (3c) 
L y = -AL z + LqL x , (3d) 
L» = AL y , (3e) 

with the spherical constraint, L 2 + Ly + t? z = 1. 

Solving the above equations set to zero yields the fixed 
points of the system. It is simple to see that there exists 
two fixed points for all parameter values, at 



±1, L z = L y = q = p = 0, 



(4) 



and for L 2 > Auj 2 there exists a further four fixed points, 
located at 




with L y = p = 0. Stability analysis of the fixed points 
shows that the original fixed points (0J are stable for 
L 2 < Au 2 , then lose their stability above this critical 
point, whilst the emergent fixed points are stable. The 
situation where a solitary fixed point becomes unstable 
and two new, stable fixed points emerge at some critical 
parameter value is called a supercritical pitchfork bifur- 
cation. The bifurcation point, L 2 < Aim 2 , corresponds to 
the critical parameter values identified in Ref. jl4| . 

The bifurcation implies that above the critical point, 
the energy is minimized by assuming a non-zero value of 
the oscillator displacement, x — ±-^L z , and the spin is 
now localized with a non-zero L z . 

In a recent paper we studied this type of bifurca- 
tion and its relationship to entanglement |16| . Lambert, 
Emary and Brandes |S| studied the entanglement in the 
more generalized system of a collection of N qubits cou- 
pled to a single oscillator. Since the qubits are all identi- 
cally coupled to the oscillator mode, they can be modeled 
as a single qudzi, meaning this system has the same clas- 
sical analogue as described in Sec. Ill Al exhibiting the 
same bifurcation. In the next section we study the qubit- 
oscillator entanglement in the ground state of Hamilto- 
nian J2J) , and finish by discussing the model of Lambert 
et al |8|. 



A. Classical analogue and bifurcations 



B. Entanglement in the ground state 



Here we clarify that the critical parameter found by 
Levine and Muthukumar [LU corresponds to a bifurca- 



In their study of the characteristics of the ground state 
entanglement between the qubit and oscillator, Levine 
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and Muthukumar [T^ focus on the determination of spe- 
cific correlation functions via functional integrals with 
the characteristics of these functions being indicative of 
entanglement. We focus solely on a quantitative study 
of the entanglement, employing the canonical measure 
of bipartite entanglement, the entropy of entanglement, 
which is the von Neumann entropy of the reduced density 
operator, p of the qubit i.e., 



S(p) = p\og 2 p. 



(0) 



To begin our study of the ground state entanglement, we 
consider the case where there is no transverse magnetic 
field applied to the qubit (i.e., A = 0). 



A = 




o o 



In the case of Hamiltonian J5J with A = the eigen- 
state problem is exactly solvable |20j . Each energy eigen- 
state is two-fold degenerate, spanned by the (orthogonal) 
states 

II) I I) (7) 

(q\€) = X n («+£)|t) =x£(9)lt) (8) 

with energies, 

L 2 

E n =ujn - —, (9) 

where Xn(s) is the n th linear harmonic oscillator wave- 
function. We see that \tp^) and \tp^} correspond to states 
localized in the left and right displaced harmonic wells, 
respectively. Note the correspondence with the fixed 
points derived earlier (|4I5[) . For the ground state, we 
have 

Xb(ff) =(£)V*«'. (10) 



FIG. 1: Entanglement in the ground state of the E®f3 system 
for different superpositions of the degenerate ground state 
(defined by Eq. lllU 'l for increasing qubit-oscillator coupling, 
L/lu. 



\i) ,|T> 




q 



From the degeneracy, a general ground state can be writ- 
ten as any superposition of the states J7J , JSJ 



l^o) = c 1 |^)+c 2 e^|Vo Ji ), 



(11) 



with c 2 + c 2 , — 1. A general density operator describing 
the ground state is thus 

p = ciXo(«) 2 l DU \ + cic 2 xo(q)xo(q) 
( e - t7 l!)(TI + e l7 IT)a I) 
+clx?(?) 2 l T)(T I- 

Tracing out the oscillator degree of freedom, the re- 
duced density operator, p s is 



1 

cic 2 e 



Q + Z7 



cic 2 e 
1 



a— 



(12) 



FIG. 2: The coupling of the oscillator displacement to the spin 
acts to split the oscillator potential in two. The ground state 
is then either localized in one of the two potential wells - states 
l^o^i l^o ) or a superposition of the two. As the coupling 
increases, the spatial separation of the two states increases. In 
turn, the overlap of the state decreases and the entanglement 
increases. The above corresponds to an equal superposition, 
which achieves the maximum entanglement. 



where a = 2(L/uj 2 ) 2 . This density operator allows the 
entropy of entanglement of the ground state to be de- 
termined as a function of c\ and the coupling L/u) (it is 
independent of the phase 7) and is shown in figure ^ 

Note there are two degenerate ground states (IV'o) an( l 
l^o 1 )); where the qubit is never entangled with the oscil- 
lator, regardless of the coupling strength. For all super- 
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FIG. 3: Entanglement in the ground state of the E®[3 system 
(with an applied transverse magnetic field) with respect to 
a = L 2 /lu 2 A for various A (with u = 1). The inset is the same 
entanglement results, but with respect to solely the coupling 
strength, L. In all cases, the entanglement becomes maximal 
as L — > oo. As the ratio of A/uj increases, the distinction 
between the "separable" and "entangled" regions identified 
my Levine and Muthukumar [l4jl , on either side of the classical 
bifurcation (a = 1), become apparent. 

positions of the two degenerate states, the entanglement 
increases, as the coupling, and hence the spatial separa- 
tion of the two states increases (see figure |2J. Maximum 
entanglement is achieved for an equal superposition. 

2. A ^ 

The addition of the Aa x term to the Hamiltonian 
means that the eigenvalue problem is no longer exactly 
solvable so the ground state must be analyzed numeri- 
cally (see appendix lA"|l . 

The Aa x term breaks the original degeneracy and 
forces the ground state to exhibit a superposition be- 
tween the 'up' and 'down' spin states, resulting in an 
entangled ground state for all A > (see figure 0). How- 
ever, a non-zero A means the oscillator potential can 
no longer be viewed as two, spatially separate harmonic 
wells each corresponding to either of the two orthogonal 
states of the spin ('up' and 'down') as shown in figure 
121 Instead, the two separated wells now each correspond 
to some superposition of the spin states. This effect re- 
sults in a decrease in the entanglement between the qubit 
and the oscillator (as evidenced in figure 0). In the limit 
of large L, the non-zero A ensures the ground state ap- 
proaches a maximally entangled, equal superposition of 
the now far spatially separated states Q) and ||SJ. 

Levine and Muthukumar focus on the entanglement 
in the ground state in the massive limit of the oscil- 
lator. More rigourously, this is defined as the limit of 
in — ► oo while keeping mtu 2 constant. Alternatively, 



this is equivalent to A/u — > oo. This is the limit of 
the quantized harmonic oscillator approaching its (con- 
tinuous) classical counterpart. In Ref. Q it is argued 
that in this massive (classical) limit of the oscillator, the 
onset of entanglement in the ground state becomes dis- 
continuous with respect to the parameter a = L 2 / Amu 2 . 
This led Levine and Muthukumar to the identification 
of "separable" and "entangled" parameter regions for the 
ground state - below the critical parameter, a c = 1, the 
ground state is separable and at a c there is a discon- 
tinuous change in the ground state, whereby it becomes 
entangled. 

In stating that the onset of entanglement becomes dis- 
continuous, it is assumed that Levine and Muthukumar 
mean that the ground state entanglement with respect 
to the parameter a becomes non-analytic at the critical 
point. This is not surprising, since in the classical limit of 
the oscillator, the ground state does not change smoothly 
with respect to a at the critical a c . 

In the classical limit, the ground state corresponds to 
the bifurcating fixed point, (L x = 1, L z = q = for a < 
1) identified in Sec. Ill Al As the oscillator behaves more 
classically, the change in the ground state with respect 
to a becomes non-analytic at a c . Due to the pitchfork 
nature of the bifurcation, the ground state transforms 
from the oscillator state localized around the single fixed 
point, to a superposition between the two emergent fixed 
points, as is passes through the bifurcation, i.e. (q) = 
for a c < 1 while (q) = — (<7 Z ) — i<7o f° r a c > 1; where 
(&z) 7^ 0. This is not the only model system where such 
a bifurcation can be used to infer an understanding of 
the entanglement properties of the ground state. 

The system considered by Lambert, Emary and Bran- 
des [i| describing the interaction of N qubits with a single 
bosonic mode Q (known as the Dicke model), undergoes 
a quantum phase transition in the N oo limit, at a 
critical value of the coupling, L — L c . Here the entan- 
glement between the ./V-qubit ensemble and the field in 
the ground state, with respect to the coupling strength, 
L was considered. It was demonstrated that the entan- 
glement obtained its maximal value corresponding to the 
critical coupling. More interestingly, the entanglement 
goes to infinity and becomes discontinuous in the N — > oo 
limit. 

The classical analogue of the Dicke model is identical 
to that defined in Sec. Ill Al with the critical coupling 
corresponding to the bifurcation in the classical analogue. 

In Ref. we demonstrated that for a system of 

coupled giant spins whose classical analogue exhibits the 
same bifurcation, the entanglement between the spins 
with respect to the coupling strength is peaked at a cou- 
pling strength corresponding to the bifurcation. In the 
limit of infinite angular momentum, the maximum en- 
tanglement goes to infinity at this critical point. 

In all three cases described above characteristics of the 
entanglement can be understood by considering the fixed- 
point bifurcation in the classical system. 

We now take the next logical step and study the ground 
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state entanglement in a system of a qubit coupled to two 
oscillators. 



III. E ® e : QUBIT WITH TWO DEGENERATE 
OSCILLATOR MODES 

The E (g) e Jahn- Teller System models the interaction 
between a doubly degenerate electronic state (E) and a 
doubly degenerate normal mode (e) 01- This is analo- 
gous to a qubit coupled to two harmonic oscillators. Fol- 
lowing the notation of Englmann, the Hamiltonian mod- 
eling this system is defined as 




H = -flLU 

2 



4 



dqj 



1 



-L(q e a e + q e a e ) (13) 



where oj is the natural frequency of the identical oscilla- 
tors and L is the vibronic coupling strength (all in units 
of H). In terms of the basis states of the qubit (or the 
electronic doublet), denoted 



14) = 



, IT) = 



o 



(14) 



the spin operators are defined as 



o-e 



-1 
1 



1 

1 



(15) 



Defining the usual oscillator mode creation and annihi- 
lation operators via 



Po 
(R 
Pi 



1 

i 



hV2 



(a + a f ) 
- ( flt " a 



hV2 



(16) 
(17) 
(18) 
(19) 



where pg = ihdg and p e = ihd t allows the Hamiltonian 
(|13fl to be written as 



H = hu (a + a + 6+6 + 1 
L 



+ 



2V2 



[(a + a t )<7fl + (6 + & t )<7 e ] (20) 



The adiabatic potential for this Hamiltonian has the 
'Mexican-hat' shape, as in figure^] Like the single oscil- 
lator case, the coupling of the qubit to the two, orthog- 
onal oscillators, results in a splitting of the no parabolic 
potential in the two spatial oscillator dimensions. 



FIG. 4: The Mexican-hat shaped potential. Vo and V\ corre- 
spond to the parabolic potentials of the individual harmonic 
oscillators, displaced from the origin by their coupling to the 
qubit. 



A. Conserved Quantity 

The total angular momentum of the system J, is the 
sum of the orbital angular momentum L (contributed by 
the harmonic oscillators) and the spin angular momen- 
tum <r (contributed by the qubit) i.e. J = L + a. 

Defining the direction q^, as that perpendicular to qg 
and q f , it is possible to show that 

J-CU L^j T 

- the total angular momentum in the nj-dircction - is a 
constant of the motion. Firstly, define 



i 
-i 



(21) 



such that <7 e ,(7 CT ,<70 correspond to the Pauli's & x ,a y ,a z 
respectively and 



Ltd = qePe - q^PB = —itl 







(22) 



Starting with the generic commutation relation relat- 
ing position, q and momentum p, [q,p] — ih, we now 
note the following commutation relations for L ro relating 
to the relevant terms in the Hamiltonian: 



L m , qe 



iq e 



-2 



2 



L^,p 



■- 2iq e 



L^,p 2 g 



= -2ip e pg. 



Together with the Pauli spin operator commutation 
relations, 



lag 



(and cyclic permutations) (23) 



6 



it is simple to see that J ro , H = 0, so J ro is a constant 

of the motion. Note that this is different to Ref. j2f| 
which claimed L ro was conserved. 



e q amh 2 (aq) q dq = —e a a-v^rErf (a) 



B. Semi-classical fixed points 

The equations of motion for the classical analogue of 
the E (g) e system are similar to those of the E ® f3 
I3ejl . except now there is an extra degree of freedom from 
the additional oscillator mode. For Hamiltonian Q13[l. 
there exists two fixed points at the origin position of the 
oscillators, with L ro = ±f, and then a ring of stable 



fixed points around the origin, satisfying L 2 
with q e = -^,L e ,q e = 
to the potential, figure 



2^tjj . Note the correspondence 



C. Ground State Ansatz 



From the work of Englman 15, 22], we now introduce 
the following ansatz for the ground state of the Hamil- 
tonian 113fl . This approximation is based on a similar 
construction to that of the eigenstates in the E ® (3 case. 



where 
A(q,cf>) 



1 



-L 2 /(2hui) 2 



(A(q^)\ i)-<S(g,0)| T» 
(24) 



B(q, . 



\2hw J 
\2fko J 

and we have adopted a polar coordinate system for the 
oscillator variables qe = qcos(4>), q e = qsm((f>). Note that 
</> commutes with q. The ground state is degenerate, and 
the orthogonal ground state to | VP) is simply its complex 
conjugate, |**), i.e., W^*) = 0. 

It was shown in [nj |23| that this ansatz gave a good 
approximation to the ground state energies of the Hamil- 
tonian In this section we shall use it to derive an 
expression for the ground state spin-oscillators entangle- 
ment. We find good agreement between this expression 
and numerical results. 

Entanglement between the spin and the two oscilla- 
tors can be calculated from the von Neumann entropy of 
the spin's reduced density matrix obtained by taking the 
partial trace over the oscillator variables. 



Ps 



\*(q,<l>))(*(q,<l>)\qd<l>dq (25) 



In calculating p s we will make much use of the following 
integrals 



f°° — 2 2 1 , 

I e 9 cosh (aq) qdq = — [ 2 + < 



^Erf(a)) 



f°° _ 2 1 2 

/ e 9 cosh(a<7) sinh(ag) q dq = -e a a^/n 
Jo 4 

where Erf(x) is the error function ranging between and 
1. We will further require the integrals 



2tt />oo 



JO 



A(q, <p)B(q, </>)* dtp qdq = ir 



The state (q,9\^)} is not normalized. Its normalization 
factor N is given by 

iV 2 = (*(g,#|tf(g,0)) 

! e -2L 2 /{2Hu>) 2 

2 

r27r poo 

x / / \A(q,<p)\ 2 + \B(q,ct>)\ 2 d<Pqdq 
Jo Jo 



-2L 2 /(2hu)) 2 



1 + //W ! ^Erff|) 
2hiu \2hLo) 



The reduced density matrix of the spin system is then 

p -2L 2 /(2hw) 2 r 2rr /-co 

Q fc- 



P 



2N 2 



(l^0)| 2 U)(l 



+A( q ^)B(q^)*\ 4)(T | 
+A(q,<f>)*B(q,<P)\ T)(| | 

+ \B(q^)\ 2 \ TXT \)dcf>qdq 
and using the above integrals this evaluates to 



1 



I iC{L/uj) 
-iC(L/uj) 1 



where C{L/uj) 



(26) 



(27) 



It can 



+ e^)^ E rf(db) 
be readily seen for large coupling L/lu — > oo we have 
C{L/w) — > and the state p s is completely mixed. The 
entanglement of formation between the spin and the os- 
cillators, given by S{p s ) takes its maximum value of 1 in 
the strong coupling limit. On the other hand for small 
coupling L/lu « f, we find C(L/oj) is also close to one 
and ps approaches a pure state and the entanglement of 
formation for the system approaches zero. 

The reduced density matrix for the orthogonal degen- 
erate ground state \ty*) is simply the adjoint of p s and 
its entanglement properties are identical. Somewhat sur- 
prisingly however, a ground state superposition of these 
two displays different entanglement properties. 
Consider an arbitrary such superposition 



I*) 



(28) 
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where c\ + c?, = 1 . Neglecting normalization for the mo- 
ment the density matrix entries for the system can be 
written in the | J.), | |) basis as 

,2 



Entanglement in Superpositions of Jahn-Teller Ground States 



Poo (<?,</>) 



= \ciA(q, 



c 2 e 



-*7 



i (ciA(q,(f)) + c 2 e' 
x ( Cl B(q,(t)) 
Pw(M,4>) = -i{ciA{q,4>) 
x (ciB(q,(f>) 
P\\(M,4>) = \ciB(q,(j)) - c 2 e n B(q,< 



B(q,4>)) 



- c 2 e 

+ c 2 e-^A(q n )) 
c 2 e^B{q, 1 Y) 



and, as before, we can calculate the reduced density ma- 
trix entries 

(■2T! 

poo (q, 4>) # dq 



Poo 



2cic 2 



2tt 







\A(q^)\ 2 dcpqdq 
Re(e- n A(q,(j)) 2 ) d(j)qdq 



(29) 



The first term we have already calculated and we find 

r 2v 



JO 

n 



Re (e'^Aiq, 0) 2 ) d(j>qdq 

L 2 /(2huj) 2 



cos(7) 1 + 



1 



Erf (JL\ 
2hw \2hw) 



Reintroducing the normalization factor into equation 1291 
gives us 



1 



Poo = 2 (! + c i c 2 cos (7) (1 +C(L/w))) 
with C(L/w) as before. A similar calculation finds 



(30) 



Poi = 2 ( c i 



4)C(L/w) 
1 



Cl c 2 sin( 7 )(l+C(L/u;)) (31) 



and since p s is a density matrix the remaining two entries 
are pf Q = (/&))*, Pn = 1 - Poo- 

The eigenvalues of ps can be written as | (l ± \/l — T\ 
where 



r = 1 - c\c\ (1 + C(LM) 2 - (c? - cl)Y C(L/lo) 2 (32) 

Interestingly we see that if we take an equal superposi- 
tion ci = c 2 = -j= then let the coupling become very 

strong L/uj — > oo these eigenvalues become j, | and the 
entanglement of formation is S(p s ) « 0.8113. The spin- 
oscillator entanglement in this ground state can never 
reach a maximum value regardless of how large the cou- 
pling term is. This is quite different from the correspond- 
ing results in the E(&(3 model, where there exists a ground 
state separable for all couplings. This will be discussed 
further in Sec. ITTTeI 




FIG. 5: The entanglement in the (ansatz) ground state of 
the E ® e model with respect to the coupling strength for 
all possible superpositions of the two degenerate states. Note 
that this is quite different from the corresponding results for 
the E ® f3 model (shown in figure . Here the ground state 
is always entangled, regardless of the superposition. 



D. Numerical analysis 

We now compare the ground state entanglement re- 
sults from the ansatz with exact numerical results. The 
Hilbert spaces of the two oscillators were truncated to 50 
basis states. Increasing the Hilbert space further had no 
effect on the results. 




FIG. 6: Comparison of the entanglement in the ground state 
between the exact results obtained by numerical diagonaliza- 
tion and the ansatz of Eq. I|24jl . 



We see that there is good agreement between the exact 
numerics and the ansatz, particularly in the small and 
large coupling limits. 
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FIG. 7: Comparison of the entanglement in the ground state 
between the numerical diagonalization and the ansatz for an 
equal superposition of the two orthogonal states. 

E. Distributed Entanglement 

The entanglement in the ground state we have con- 
sidered so far is that between the qubit and the pair of 
oscillators. Since the oscillators, and the couplings to the 
qubit, are identical, the entanglement to the qubit is dis- 
tributed equally between the two oscillators. However, 
it is possible to consider quantitatively how the entan- 
glement between the qubit and the oscillators is shared 
between the two polar degrees of freedom - radial and an- 
gular coordinates. Such entanglement involving (orthog- 
onal) internal degrees of freedom, as opposed to physical 
partitions of the system have been considered, for exam- 
ple, in the context of trapped ions p^ |. where spin and 
orbital degrees of freedom of a single ion are entangled. 

In the limit of large coupling, L/lo ^> 1, the ground 
state ansatz can be expressed as 

(q, m « H<l) (10) - + e» (|0) + (33) 

where T{q) = e (-£ 2 /(2M 2 e ( 9 L/2M- g 2 /2/(^2jv). Thc 
radial coordinate, q, is separable, hence the qubit is en- 
tangled solely with the angular degree of freedom, 4> of 
the two oscillators. Outside of this parameter range how- 
ever, the radial coordinate is not separable, meaning the 
qubit is entangled with both degrees of freedom. 

To quantify this distribution of entanglement for thc 
ground state H24jl. it is possible to determine the entan- 
glement solely between the angular degree of freedom and 
the qubit. We begin by identifying the states 

U m {<l>) = -L e ±^ (34) 

as eigenstates of L ro = ihd^, with eigenvalue m. In the 
ground state, \ip), only the m — 0, 1 states are present. 
So the angular degree of freedom, 0, is constrained to a 



0.09 




0.5 1 1.5 2 2.5 3 3.5 4 

L 



FIG. 8: The difference between the qubit-oscillators entangle- 
ment and the qubit-angular degree of freedom entanglement 
(in terms of the von Neumann entropy) as a function of L/u. 
Hence, almost all the entanglement is between the qubit and 
the angular degree of freedom. 

two-dimensional subspace of its total Hilbert space. Let- 
ting 1 0) = Uq{4>) and |1)^ = Ui(4>), we may view the 
angular degree of freedom in the ground state as itself a 
qubit, reducing the problem of thc entanglement between 
the (spin) qubit and </>, to the well-known situation of a 
pair of qubits. Rewriting the state of the(spin) qubit in 
the basis |+> = (| |) + i\ T» /y/2 R = (| |> - i\ ?)) /y/2 
(which are the eigenstates of cr ro ), the ground state in the 
limit of large coupling, Eq. I|33[) becomes 

(g|*)«V^te)(|->|o>, + |+)|i>,) (35) 

Clearly, the spin-qubit and the '<^'-qubit are in a max- 
imally entangled Bell state, completely separable from 
the radial coordinate. 

The concurrence p3 . l25j | is a good measure of the two- 
qubit mixed-state entanglement, which we can use to 
quantify the entanglement between the spin-qubit and 
the <^>-qubit. 

The concurrence, C, between a pair of qubits, A and 
B is defined using the "spin-flipped" density matrix 

Pab = [fry ® a-y) p* AB (a y ® a y ) (36) 

where the asterisk is the complex conjugation in the stan- 
dard basis. If the square roots of the eigenvalues of the 
product pab Pab in decreasing order are L±, L2, £3, £4, 
then the concurrence of the density matrix pab is 

C = min{0,Li- L 2 - L 3 - U} . (37) 

The concurrence is related to the von Neumann en- 
tropy via the tangle, r = C 2 , by 

S -(i±£Hi) 



9 



where H(x) = — x\og 2 (x) — (1 — x)log 2 (l — x) is the 
Shannon entropy. 

Using the above, it is possible to calculate the entropy 
of entanglement between the qubit and the angular de- 
gree of freedom, and compare it to the total entangle- 
ment between the qubit and the two oscillators. Figure 
[H] shows the difference, AS, between these two entangle- 
ments. AS asymptotes to zero, such that in the strong 
coupling regime, the qubit becomes disentangled from 
the radial degree of freedom and is solely entangled with 
the angular degree of freedom, as predicted by Ea. (|551) . 
Furthermore, this entanglement is maximal. Note that 
AS is relatively small implying that qubit-oscillators en- 
tanglement is concentrated between the angular degree 
of freedom and the qubit. 

The two orthogonal degenerate ground states, and 
it's complex conjugate, from the ansatz Eq. (|24|l . 
are the two sole ground states where the angular degree 
of freedom can be treated analogous to a qubit. In any 
superposition of these two states, the states of the angu- 
lar degree of freedom is in the subspace spanned by the 
states lA m {4>) with m — 0, ±1 - now a three-level system, 
or qutrit. 

As shown in figure^ for any superposition, the ground 
state entanglement does not asymptote to the maximal 
value, but however, there is no ground state superposi- 
tion that has zero entanglement, as in the single oscillator 
(E <g> (3) case. In the large coupling limit, the radial de- 
gree of freedom still becomes separable, such that the 
entanglement is concentrated between the qubit and the 
angular degree of freedom for all ground possible states. 
The observation that in all superpositions the angular de- 
gree of freedom is viewed as a qutrit rather than a qubit 
could explain why the entanglement is never zero - as 
seen in the single oscillator case - nor maximal. 



F. Addition of transverse magnetic field 

For completeness, we now consider the effect of ap- 
plying a transverse magnetic field to the qubit, in the 
direction perpendicular to the oscillator displacements. 
The Hamiltonian thus becomes 



H = Ac 



1 / 2 



Qe 



Pi 



1 



L {qecre 



<lsOe) 

where A is the strength of the magnetic field. 

With respect to the fixed point structure in the clas- 
sical analogue, the addition of the A term has the effect 
of destroying the stable ring of fixed points, leaving four 
stable points, at 



L P = ±\ l- 



16w 2 A 2 



L 1 



4cjA 



and 



16w 2 A 2 



L 4 



4wA 



This implies that again a pitchfork bifurcation is present, 
with the critical coupling, L 2 — 16a-> 2 A 2 . This should 
again manifest itself in the large A limit of the entangle- 
ment in the ground state, as in the E <g> (3 model. 

Moving to the quantum regime, the ground state is 
no longer degenerate, and the presence of the transverse 
field forces the ground state into a maximally entangled 
state in the large coupling limit. 
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FIG. 9: Entanglement in the ground state of the E($s model 
in the presence of a transverse magnetic field. Note the simi- 
larities with the E <g) j3 model results in figure [3] 

From figure^ it's clear that the bifurcation in the E(g>e 
model plays a similar role as that in the E ® (3 model, in 
the large A/ui limit (the classical limit of the oscillator). 
The gradient of the entropy of entanglement curve with 
respect to L becomes steeper around the critical point, 
and we see the division into the "separable" and "entan- 
gled" parameter regions. 



IV. CONCLUSION 

We have studied the entanglement in the ground states 
of the E (g) (3 and the E ® e Jahn- Teller systems, which 
model a single qubit coupled to one and two harmonic 
oscillators, respectively. 

In the single oscillator case, we have considered the re- 
sults of both Levine and Muthukumar 0] and Lambert 
et al. ||. In both cases, we have argued that the entan- 
glement characteristics of the ground state can be under- 
stood by considering the bifurcation of the fixed points 
in the classical counterpart. In the two extremes consid- 
ered in and ||| the classical limit becomes relevant - 
either of the oscillator, or the entire system, respectively. 
Again, as shown in previous work |lq | , the nature of the 
bifurcation (the pitchfork structure) is crucial - a single 
fixed point becomes two, leading to superposition states 
in the quantum regime. 
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In the E ® e model, we found that the ground state 
entanglement between the qubit and the oscillators dif- 
fered from that of the single oscillator model, insofar as 
for no superposition of the orthogonal ground states was 
there zero entanglement. Furthermore, how this entan- 
glement is shared between the two degrees of freedom of 
the double oscillator subsystem was considered. It was 
found that the entanglement between the qubit and the 
two oscillators is concentrated between the qubit and the 
angular coordinate, with the radial coordinate becoming 
completely separable in the large coupling limit. This 
correlation between the angular degree of freedom and 
the qubit states is not surprising given the radial sym- 
metry of the potential created by the qubit-oscillators 
coupling. 

The Hamiltonian of the E<g>e model in Eq. (|2C)|> can be 
physically realized using two vibrational degrees of free- 
dom of a single trapped ion .26|. The required coupling 
is achieved using external laser pulses to couple different 
components of the atomic polarization vector, a to each 
of the vibrational modes. 

In Ref. 01 1 Sjovist used the E®e Jahn- Teller system 
as a model for electron nuclear interaction. While the 
entanglement in higher energy eigenstates was considered 
in that article, our results for the ground state in the large 
coupling limit coincide. Our results will hopefully shed 
more light on the characteristics of this electron-nuclear 
entanglement in molecular ground states. 

One of the most intriguing results of this paper is that 
the E ® e Jahn- Teller model always has an entangled 
ground state and when we take the semi-classical limit 
(L/uj — > oo) the entanglement between the qubit and the 
oscillators approaches its maximal value. In contrast, for 
the E (g> (3 model there are two degenerate ground states 
for which there is no entanglement. It appears that this 
difference is due to the presence of the angular degree 
of freedom for the oscillators. We conjecture that the 
entanglement is intimately connected with the geometric 
(Berry's) phase associated with cyclic adiabatic varia- 
tions of the angular co-ordinate of the classical limit of 
this model p7| . 

The above raises an important question as to whether 
our results are a manifestation of a very general phenom- 
ena connecting entanglement and geometric phases. In 
the hope of stimulating further work we offer the follow- 
ing conjecture. 

Conjecture: Let H(S) be a Hamiltonian which de- 
pends on some parameter S and acts on a bipartite 
Hilbert space V = Vq ® Vq of 'finite dimension. Suppose 
that in some limit S — > S c i, the Hamiltonian becomes 
H(C) which acts on the Hilbert space Vq where C de- 
notes a finite dimensional parameter. Suppose also that 
there is geometric phase associated with cyclic adiabatic 
variations of C. Then for all possible ground states of 
H(S) there is always entanglement between Vq and Vc- 
Furthermore, the entanglement approaches its maximum 
possible value as S — > S c i . 

This conjecture should first be tested for the T <g> H 



Jahn- Teller model which describes three-fold degenerate 
electronic levels coupled to a five-fold degenerate phonon 
and which is relevant to fullerene (Ceo) molecules psj . 
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APPENDIX A: NUMERICAL BASIS - 
DISPLACED FOCK STATES 

Numerical analysis of a system within an infinitely 
dimensional space often implies some truncation of the 
Hilbert space for calculations. For the E ® (3 system, to 
reduce the potential numerical error from this truncation, 
rather than choosing the set of Fock states as the basis for 
the Hilbert space of the oscillator, we use the displaced 
Fock states (JU,© corresponding to the eigenstates for 
A = 0. In this basis, the Hamiltonian is diagonal for 
A = 0, with entries given by the energy eigenvalues 
For, A ^ 0, we must calculate the off-diagonal elements. 
Since this set of basis states is not orthogonal, to deter- 
mine the off-diagonal elements of the Hamiltonian matrix 
we make use of the following expressions: 

(V>mW> = (^ml$?> = 5 mn, 

/ Xm(q)x%(q)dq = s mn , 
I Xm(q)x%(q)dq = s mn , 

fxtk(q)x%(q)dq=(m\D(2L/oj 2 )\n), 

J X*(q)Xn(l)dq = (m\D(-2L/u 2 )\n), 

where D(a) is the displacement operator. From Caves et 
al. 123 we have 



(m\D(/3)\n) 



n\ am—n r ( m_ n ) 

HZ 

e 2 



7rf(/3*) r 



' L 



(n — 7n) 



(W\ 2 



if m > n 



if n > m 



(_l)m+™ e -2- 

where L^.(u) are the generalized Laguerre polynomials. 
Now, considering N + 1 oscillator modes, the non-zero 
off-diagonal elements of Hamiltonian (J2J matrix, TC are 
given by 

H(m+l,n + N + 2) = (^\H\^) 

= A(xM(T |*x|i) 

= A(m\D"< {-L/uj 2 ) D {L/uj 2 ) \n) 

= A(m\D (2L/u 2 ) \n) 
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and similarly, 

H(m + N + 2, n + 1) = A(m\D (-2L/lu 2 ) \n), 



[1] H.L. Haselgrove, M.A. Nielsen and T.J. Osborne (2003), 

quant-ph/0303022. 
[2] J.I. Latorre, E. Rico and G. Vidal (2003), quant- 

ph/0304098. 

[3] T.A. Costi and R.H. McKenzie, Phys. Rev. A 68, 034301 
(2003). 

[4] A.P. Hines, R.H. McKenzie and G.J. Milburn, Phys. Rev. 

A 67, 013609 (2003). 
[5] A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 

416, 608 (2002). 
[6] T.J. Osborne and M.A. Nielsen, Phys. Rev. A 66, 032110 

(2002). 

[7] D. Gunlycke, S. Bose, V.M. Kendon and V. Vedral, Phys. 

Rev. A 64, 042302 (2001). 
[8] N. Lambert, C. Emary and T. Brandes (2003), quant- 

ph/0309027. 

[9] H.L. Haselgrove, M.A. Nielsen and T.J. Osborne (2003), 
quant-ph/0308083. 
[10] F. Verstraete, M.A Martn-Delgado and J.I. Cirac (2003), 
to be published in Physical Review Letters, quant- 
ph/0311087. 

[11] G. Vidal, J.I. Latorre, E. Rico and A. Kitaev, Phys. Rev. 

Lett 90, 227902 (2003). 
[12] A.J. Leggett, S. Chakravarty, A.T. Dorsey, M.P.A. 

Fisher, A. Gerg and W. Zwerger, Rev. Mod. Phys. 59, 1 

(1987). 

[13] U. Weiss, Quantum Dissipative systems (World Scientific, 
1999), 2nd ed. 

[14] G. Levine and V.N. Muthukumar (2003), cond- 
mat/0301463. 



both of which can be evaluated using the expressions 
above. 



[15] R. Englman, The John-Teller Effect in Molecules and 
Crystals (John Wiley & Sons Ltd., 1972). 

[16] A.P. Hines, G.J. Milburn and R.H. McKenzie (2003), 
quant-ph/0308165. 

[17] E. Sjoqvist, Int. J. Quantum Chem. 77, 526 (2000). 

[18] C. Emary and T. Brandes (2003), cond-mat/0301273. 

[19] P. Glendinning, Stability, instability and chaos: an intro- 
duction to the theory of nonlinear differential equations 
(Cambridge University Press, 1994). 

[20] G.D. Mahan, Many-particle physics (Plenum Press, New 
York, 1990), 2nd ed. 

[21] H. Yamasaki, Y. Natsume, A. Terai and H Nakamura 
(2003), cond-mat/0304484. 

[22] R. Englman, Phys. Letters 2, 228 (1962). 

[23] C. Monroe, D.M. Meekhof, B.E. King, W.M. Itano and 
D.J. Wineland, Phys. Rev. Lett. 75, 4714 (1995). 

[24] S. Hill and W.K. Wooters, Phys. Rev. Lett. 108, 171 
(1997). 

[25] W.K. Wooters, Q. Inf. Comp. 1, 27 (2001). 

[26] J. Steinbach, J. Twamley, and P. L. Knight, Phys. Rev. 
A 56, 4815 (1997). 

[27] C. A. Mead, Rev. Mod. Phys. 64, 51 (1992). 

[28] J.E. Han and O. Gunnarsson, Physica B 292, 196 (2000). 

[29] CM. Caves, K.S. Thorne, R.W.P. Drever, V.D. Sandberg 
and M. Zimmermann, Rev. Mod. Phys. 52, 342 (1980). 

[30] When there is no coupling of electronic and nuclear 
states, the problem of molecular structure is greatly sim- 
plified by the so-called Born-Oppenheimer approxima- 
tion 



